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ABSTRACT 

Studies of absorption spectra of high-z QSOs have revealed that the intergalactic medium at z ~ 2 — 3 
is enriched to ~ 10~ 3 — 10 -2 Zq for gas densities more than a few times the mean cosmic density, but 
have not yet produced an accurate metallicity estimate, nor constrained variations in the metallicity with 
density, redshift, or spatial location. This paper discusses the "pixel optical depth" (POD) method of 
QSO spectrum analysis, using realistic simulated spectra from cosmological simulations. In this method, 
absorption in Lya is compared to corresponding metal absorption on a pixel-by-pixel basis, yielding 
for each analyzed spectrum a single statistical correlation encoding metal enrichment information. Our 
simulations allow testing and optimization of each step of the technique's implementation. Tests show 
that previous studies have probably been limited by C IV self-contamination and O VI contamination 
by H I lines; we have developed and tested an effective method of correcting for both contaminants. We 
summarize these and other findings, and provide a useful recipe for the POD technique's application 
to observed spectra. Our tests reveal that the POD technique applied to spectra of presently available 
quality is effective in recovering useful metallicity information even in underdense gas. We present an 
extension of the POD technique to directly recover the intergalactic metallicity as a function of gas 
density. For a given ionizing background, both the oxygen and carbon abundance can be measured 
with errors of at most a factor of a few over at least an order of magnitude in density, using a single 
high-quality spectrum. 

Subject headings: cosmology: theory — intergalactic medium — quasars: absorption lines 



1. INTRODUCTION 

Studies of absorption lines in QSO spectra have es- 
tablished that at high redshift [z ;> 2) the intergalac- 
tic medium (IGM) is enriched by metals to metallicity 
10 -3 <; ZJ Zq <; 10 -2 at densities above a few times the 
cosmic mean (e.g., Tytler et al. 1995; Cowie et al. 1995; 
Songaila & Cowie 1996; Rauch, Haehnelt, & Steinmetz 
1997; Dave et al. 1998). Since stars are thought to form in 
appreciable quantities only inside much denser mass con- 
densations, the observed enrichment indicates that galax- 
ies, which form from the IGM, also feed back some of their 
nuclcosynthetic products in a process which is currently 
poorly understood but which is likely to be a crucial in- 
gredient of galaxy formation and evolution. 

Because galaxies form in relatively high-density regions, 
metals in low-density regions must have traveled a signif- 
icant distance. This puts strong constraints on the in- 
tergalactic (IG) enrichment mechanism, therefore deter- 
mining the IG metallicity at the lowest possible densities 
is crucial for understanding the feedback of metals from 
galaxies to the IGM (see, e.g., Aguirre et al. 2001; Madau, 
Ferrara & Rees 2001). 

Previous studies, performed primarily with the Keck 
HIRES instrument, have used measured metal line col- 
umn densities to show that, at z « 3, most absorbers 
with logA^(HI) > 14.5 (corresponding to gas overdensi- 
ties 6 = p/(p) ^5 — 10 [Schaye 2001]) have associated 
carbon lines (e.g., Tytler at al. 1995; Cowie et al. 1995; 



Songaila & Cowie 1996; Ellison et al. 2000). These studies 
have not, however, made clear how this metallicity changes 
with spatial location, redshift, or gas density. Studies em- 
ploying line-fitting could in principle treat the first two is- 
sues. But pushing to lower gas densities where individual 
metal lines are undetectable requires a statistical analysis; 
at high redshift this is also required for metals absorp- 
tion falling blue-wards of Lya, where the metals lines are 
difficult to disentangle from H I absorption. 

Ellison et al. (1999, 2000; see also Cowie & Songaila 
1998) have compared two ways of doing this: first by stack- 
ing the spectra of the metal-line regions corresponding to 
a set of low column density Lya lines, or second by com- 
paring the Lya optical depth in each pixel with the optical 
depth in the pixel at the corresponding metal-line wave- 
length. The latter "pixel optical depth" (POD) method, 
which was devised and first used by Cowie & Songaila 
(1998), has several advantages: it is fast, it is objective, 
it preserves more of the spectrum's information, it can be 
applied to heavily contaminated regions, and it appears to 
be more sensitive to metals in low-density gas, even in un- 
contaminated regions. Its shortcoming, however, has been 
that it is not immediately clear how to interpret the results 
of the POD method, and it is unclear how to convert the 
set of PODs into a metallicity Z at a given overdensity 5, 
especially when contamination is important, as for O VI 
(see Schaye et al. 2000a). 

The purpose of the present study is to present several 
important improvements to and extensions of the POD 
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technique, to provide a complete and useful explanation 
of the method, and to test it using realistic spectra gen- 
erated from hydrodynamic cosmological simulations. The 
method as described here will be applied to a sample of ob- 
served QSO spectra in a future paper. The present paper 
should also serve as a useful reference for prospective users 
of the POD technique. The testing we perform allows us to 
estimate the accuracy with which H I and metal line PODs 
can be recovered (and to refine the method for doing so), 
and to draw the connection between the optical depths 
and Z{8) for the gas. By comparing recovered quantities 
directly to the true (simulated) ones, we demonstrate that 
the POD method does work, and that we can interpret 
directly the meaning of features in the POD statistic. 

Section || contains the basics of the method, pointing 
forward to figures in the main text, for illustration. This 
section is meant for readers not interested in all the details, 
or as a general overview for readers unfamiliar with pixel 
statistics. The full details of the method, and of the tests 
we performed, can be found in the following sections. In 
§ [| we describe the simulations used in this study and how 
we generate spectra from them. In § ^ we give a careful 
account of the recovery of metal and H I PODs from QSO 
spectra, testing the individual steps by comparing recov- 
ered and "true" optical depths; this section can be used 
as a reference for the POD method of spectrum analysis, 
and is summarized in the Appendix. Section ^ discusses 
and t ests the inte rpretation of the recovered PODs: Sec- 
tions 5.1 and 5.2 relate the H I and metal PODs to the 



density of the absorbing gas, and § 5.3 shows the rela- 



tion between metal PODs and the assumed gas metallic- 
ity, given various levels of noise and other uncertainties. 
In § H we formulate and test a procedure whereby Z(S) 
can be recovered directly and perform simple tests of this 
procedure. Finally, we discuss our results and conclude in 
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2. OVERVIEW OF THE METHOD 



The aim of the method is to statistically obtain the phys- 
ical properties of the gas that produces the absorption in 
a given redshift range. The first step is to infer, in the 
presence of noise, contamination, saturation, etc., for each 
spectrum pixel the best possible estimate of the optical 
depth TLy Q due to Lya aborption by gas at some redshift, 
as well as the optical depths Tz for various metal transi- 
tions in gas at that same redshift. The second step con- 
verts the inferred optical depths to physical gas densities, 
using the tight relation that exists between them. This 
yields density ratios of H I to various metal species such 
as C IV and O VI. Ultimately, however, we aim to measure 
the metallicity as a function of overdensity. So thirdly, we 
use ionization corrections (as determined by the photoion- 
ization package CLOUDY), for gas at given density and 
temperature illuminated by a chosen ionizing background. 

Simulated spectra. We use simulated spectra generated 
from hydrodynamical simulations to test and illustrate the 
method. We patch together the physical state of the gas 
along many uncorrelated sight-lines through the simula- 
tion box, in order to obtain one long sight-line of redshift 
extent z ~ 1.5 — z qso , where z qso is the assumed emis- 
sion redshift of the quasar. (We illustrate the method for 
z qso = 2.5 and 3.5). Metals are distributed by hand, for ex- 
ample by assuming a uniform metallicity of 1% solar. We 



then generate simulated spectra (using our chosen ionizing 
background to compute ionization balances), taking into 
account all hydrogen Lyman transitions and absorption by 
important metal transitions. The wavelength extent and 
signal-to-noise ratio of these simulated spectra are chosen 
to be similar to the very high quality echelle spectra for 
bright quasars, as obtained by HIRES on Keck, or UVES 
on the VLT. Example spectra are shown in Figures 0-0. 

Recovering the optical depth. The primary issue in re- 
covering the optical depth r from the flux F = exp(— r) 
is that information is lost when the optical depth is ei- 
ther too high, too low, or contaminated by other transi- 
tions. For hydrogen, we use the higher order transitions 
to recover the optical depth in saturated pixels. However, 
these higher order pixels are (potentially) contaminated 
by other hydrogen transitions from lower redshifts. (Note 
that our simulated spectra have a realistic level of this type 
of contamination, because they are generated from realis- 
tically long sight-lines through simulations that reproduce 
the observed evolution of the Lya forest.) For hydrogen, 
we therefore consider only those pixels of the Lyman series 
which are neither too black, nor too close to the contin- 
uum. The optical depth of that redshift pixel is then the 
minimum optical depth of all the transitions that qualify, 
thereby minimizing contamination both by Lyman series 
lines from lower redshift, by from other transitions. Fig- 
ure [j] shows how well the recovered optical depth matches 
the true one. 

For metal transitions occuring in the Lya forest, cor- 
recting for contamination is of course crucial. For O VI 
and N V, we use the Lya optical depth determined ear- 
lier to correct for contamination by higher-order H I lines. 
We further reduce contamination by using the minimum 
optical depth for each of the components of the doublet. 
For C IV and Si IV, it is possible - and very important 
- to correct for self-contamination (i.e., where the sec- 
ond component of the doublet happens to fall on top of 
the first component of another system), and we find that 
an iterative correction works well (see Figure (sj). On a 
pixel- by-pixel basis the recovery is good but imperfect, es- 
pecially at low optical depth; recovered O VI and N V 
PODS are also systematically high, due to contamination 
(see Fig. ||) . To obtain a statistical signal we bin the recov- 
ered metal PODs according to the corresponsing hydrogen 
optical depth, and find the median metal POD in each bin, 
as illustrated, for example, in Figures |^ and 0^ 

These binned, recovered PODs can be further improved 
using the simulations. By employing simulated spectra 
with the same properties as the given spectrum to com- 
pute the relation between "true" and "recovered" PODs, 
we can "invert" the recovered PODs (which are affected 
by contamination, noise, etc.) to obtain an estimate of 
the "true" PODs in the given spectrum. This process is 
illustrated in Fig. |l3| . 

From optical depth to density. In the simulations (and 
presumably in reality), there is a tight relation between 
optical depth and density, both for hydrogen (see figure |]) 
and for metals. These relations allow us to invert the opti- 
cal depths and obtain the underlying (redshift space) den- 
sities. 

Ionization corrections. To convert absorption densities 
to physical densities, we require ionization corrections, 
which depend on the shape and amplitude of the ionizing 
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background, as well as on the density and temperature. 
Fortunately, the latter two quantities are related in the 
simulations, and our simulations produce a temperature- 
density relation that is consistent with the relation mea- 
sured by Schaye et al. (2000b) from the widths of Lya lines 
as a function of their column density. Therefore, assum- 
ing an ionizing background (normalized to match the mean 
observed absorption by the IGM), we can finally deduce 
true physical densities of various atomic species, obtain- 
ing for example the oxygen abundance as a function of the 
hydrogen (over) density. 

By applying this method to simulated spectra, we can 
estimate the importance of various physical effects and 
modeling uncertainties. For example, we investigate the 
effects of uncertainties in the temperature-density and 
density-column density relations, continuum fitting errors, 
partial wavelength coverage and noise. We demonstrate 
that despite these uncertainties, the algorithm is very pow- 
erful. In Figs. [[3] and [l4| we further show that we can ac- 
curately recover the metallicity of gas down to (at least) 
the cosmic mean gas density, given an accurate estimate 
of the ionizing bakground. 

We conclude that relative to other techniques of analyz- 
ing absorption spectra, there are several advantages to the 
POD method. First, because it is not based on line fitting, 
the analysis is extremely fast, objective, and its implemen- 
tation is straightforward. Second, it is accurate and robust 
because it uses the information of all transitions for each 
redshift pixel. Third, absorption can be detected statisti- 
cally, enabling us to probe lower optical depths, which is 
crucial for studying low levels of metal contamination. We 
now proceed to describe and test the various steps in more 
detail. 

3. GENERATION OF SPECTRA FROM SIMULATIONS 

To test the optical depth technique we require realis- 
tic absorption spectra with wavelength range, resolution, 
pixel size, and noise properties similar to the observed 
spectra. Besides absorption by Lya and the metal transi- 
tions of interest, the spectra must contain absorption from 
contaminants such as higher-order Lyman lines and addi- 
tional metal lines. Before describing the procedure used 
to produce such spectra from a cosmological simulation, 
we will briefly summarize the properties of the simulation 
itself. 

In this study we employ a hydrodynamical simulation 
of the currently popular flat, scale invariant, vacuum en- 
ergy dominated cold dark matter model. The simulation 
was performed with a modified version of HYDRA (Couch- 
man, Thomas, & Pearce 1995) and uses smooth particle 
hydrodynamics (SPH). The evolution of a periodic, cu- 
bic region of the universe of comoving size 12 h^ 1 Mpc 
was followed to redshift z — 1.5. The model has a to- 
tal matter density O m = 0.3, vacuum energy density 
fl\ = 0.7, baryon density fi&/i 2 = 0.019, Hubble constant 
Ho = 65 kms" 1 Mpc -1 , and the amplitude of the initial 
power spectrum is normalized to erg = 0.9. The simulation 
uses 256 3 gas particles and 256 3 cold dark matter particles, 
yielding particle masses of 2.0 x 10 6 M Q and 1.1 x 10 7 M 
respectively. This mass resolution has been demonstrated 
to be sufficient to resolve the Lya forest spectra (Therms 
et al. 1998; Bryan et al. 1999; Schaye et al. 2000b). Gas 
particles are converted to collisionless star particles if they 



satisfy the following three conditions: (1) the density ex- 
ceeds 80 times the mean baryon density; (2) the density 
p/niH exceeds 10~ 2 cm~ 3 ; (3) the temperature is less than 
2 x 10 4 K. Feedback from star formation is not included. 
The IGM is assumed to be of primordial composition with 
a helium abundance of 0.24 by mass (metals are added only 
during the generation of spectra from the completed simu- 
lation and do therefore not contribute to the cooling of the 
gas). The gas is photoionized and photoheated by a model 
of the UV-background, designed to match the evolution of 
the mean Lya absorption and the temperature-density re- 
lation measured by Schaye et al. (2000b). 

The final spectra have characteristics typical of spec- 
tra taken with HIRES/KECK or UVES/VLT. The spectra 
span a wavelength range 3000-7500 A and include absorp- 
tion by gas from redshift z — 1.5 to the redshift of the QSO 
(zqso), which is set to either 2.5 or 3.5. Absorption from 
Lyl (Lya,A1216) to Ly31, C III (A977), C IV (AA1548, 
1551), N V (AA1239, 1243), O VI (AA1032, 1038), Si III 
(A1207), Si IV (A1394, 1403), Fe II (A 1145, 1608, 1063, 
1097, 1261, 1122, 1082, 1143, 1125) is included. These 
transitions account for nearly all absorption by z — 1.5-3.5 
absorbers in the 3000-7500 A window. Other commonly 
observed absorption lines, mainly from neutral and singly 
ionized species, occur in gas of densities higher than the 
maximum gas density in the simulation. However, with 
the exception of damped Lya absorption, the observed 
spectral filling factor of this gas is small enough that the 
effect of these absorbers is negligible for our purposes. The 
simulated spectra include absorption from absorbers with 
z > 1.5 only. However, contamination by absorbers with 
redshift z < 1.5 is minimal because nearly all commonly 
observed absorption lines have rest wavelengths in the far 
UV. 

At many hundreds of snapshot times, the position, den- 
sity, temperature and velocity of those gas particles that 
affect the physical state of the gas along 6 randomly chosen 
sight-lines through the simulation box (i.e., those particles 
whose SPH smoothing kernel intersects the sight-line) are 
saved. The time between outputs is similar to the time it 
takes light to cross the simulation box. A long, continuous 
sight-line spanning z = 1.5 to z qso is created by connecting 
short sight-lines, choosing at random from those sight-lines 
with redshifts close to the desired redshift. Before connect- 
ing a short sight-line it is cycled periodically so that the 
point at which the inferred Lya transmission is maximum 
is at the ends. Since the simulation box is sufficiently 
large for each short sight-line to have a maximum in the 
Lya transmission close to 1 (i.e., no absorption), there are 
generally no strong discontinuities in the long sight line. 
After a short sight-line has been inserted in the continuous 
sight-line, the gas density along the sight-line is scaled to 
account for the Hubble expansion during the time it takes 
light to cross the simulation box (this scaling is not very 
important because the simulation box is small). Correla- 
tions in physical quantities on scales greater than the box 
size are of course not modeled correctly. However, for the 
purpose of studying metal absorption all that matters is 
that the contamination by absorbers at various redshifts 
is modeled realistically. 

The ionization balance of each gas particle is com- 
puted using the publicly available photoionization package 
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CLOUDY 1 (version 94; see Ferland 2000 for details), as- 
suming the gas to be optically thin. The gas is illuminated 
by the redshift-dependent model of the UV/X-ray back- 
ground of Haardt & Madau (1996), updated in Haardt & 
Madau (2001) 2 , which includes contributions from QSOs 
and galaxies. The amplitude of the background radiation 
is scaled so that the mean Lya absorption matches the 
measurements of Schaye et al. (2000b) 3 . The UV/X-ray 
background used to compute the spectra is different from 
the model of the background radiation that was used in the 
hydrodynamical simulation. The exact UV/X-ray back- 
ground used in the simulation is, however, unimportant as 
long as the resulting thermal evolution of the IGM roughly 
matches the observations, which is indeed the case. The 
gas fraction in each ionization state (e.g., C IV/C) depends 
on the density, temperature, and redshift. 

For each short sight-line the absorption spectra of the 
various transitions, i.e., the optical depths as a function 
of redshift, Ti(z), are computed using the formalism de- 
scribed in Theuns et al. (1998). For each transition i with 
rest wavelength A^ , the spectrum n (z) is then added to the 
final, continuous spectrum t(A) using A = A$(l + z). The 
final, continuous spectrum r(A) is converted into a nor- 
malized, flux spectrum using F = exp(— r). It is then 
processed in three steps to give it characteristics simi- 
lar to observed spectra taken with HIRES/KECK and 
UVES/VLT. First, it is convolved with a Gaussian with 
FWHM of 6.6 kms -1 to mimic instrumental broadening. 
Second, it is rcsamplcd onto 0.04 A pixels. Third, noise 
is added. For simplicity we have assumed the noise to be 
Gaussian and to consist of two components, one with flux- 
independent rms amplitude (1/3)(S/N)~ 1 and one with 
amplitude (2/3)F(X)(S/N)- 1 , where S/N is the signal- 
to-noise ratio in the continuum. For most tests we quote 
results with both S/N = 25 and S/N = 100. Of course, 
when comparing with real data one should impose noise 
that is statistically equivalent, both in its dependence on 
wavelength and flux, to the noise in the observations. 

To facilitate a physical interpretation of the results pre- 
sented in this work, we have also computed the gas density 
along the continuous sight-line. Because of redshift space 
distortions, multiple gas elements can contribute to the 
optical depth of a given pixel. We therefore follow Schaye 
et al. (1999) and use optical depth weighted physical quan- 
tities. The density of a pixel in velocity space is defined as 
the sum, weighted by the contribution to its optical depth, 
of the density of all the gas elements that contribute to the 
absorption in that pixel by a given transition. Note that 
because different elements have different ionization correc- 
tions, the density weighted by, say, the C IV optical depth 
will generally differ somewhat from the density weighted 
by the H I optical depth, even if the mctallicity is uniform. 
This is not an artifact of the density definition; there is no 
reason why the C IV and H I absorption at a given redshift 
should arise in exactly the same gas. 



4. THE OPTICAL DEPTH TECHNIQUE 
1 http : / /www . pa . uky . edu/ ^gary/ cloudy/ 

2 The data and a description of the input parameters can be found e 

3 Note that the standard practice of rescaling the spectrum (i.e., mul 
the background radiation, could lead to large errors if absorption by m 



We now describe the POD analysis method we have 
tested using the simulated spectra generated as described 
in § |3|; real spectra can be analyzed in the same man- 
ner after continuum fitting and rescaling. The basic tech- 
nique was pioneered by Cowie & Songaila (1998; see also 
Songaila 1998) and improved by Ellison et al. (2000) and 
Schaye et al. (2000a). In the POD method, the metal- 
line optical depth in each pixel is compared to the opti- 
cal depth in the pixel corresponding to H I absorption by 
gas at the same redshift. The array of pixel pairs is then 
binned in H I POD, and the median metal POD in each 
bin computed. For a uniform density ratio of the metal 
ion to H I, and with perfect POD recovery, the result- 
ing relation would be linear. In general, the detection of 
the metal species will manifest as a correlation between 
metal POD and H I PODs, and the change in the slope of 
this correlation indicates changes in the ratio of densities 
(resulting either from changes in metallicity or ionization 
balance with density). The utility of the POD technique 
relies upon an understanding of how a given physical sit- 
uation maps into the recovered POD relation. Here we 
review the details of the method, and while doing so, com- 
pare 'recovered' (after noise, blending, etc.) quantities to 
'true' (i.e., simulated) quantities where possible. We give 
a concise summary of our recommended implementation 
of the method in the Appendix. 

Given a QSO redshift z qso , we identify the redshift 
zp = (1 + 2qso)(ALy/3/ALya) — 1 such that absorption by 
gas in the redshift interval \zp^ Zq SO ] is uncontaminatcd by 
Ly/3 absorption. We will analyze pixels in the redshift 
range [z min , z max ] with zp < z min < z max < z qso . For a 
transition with rest wavelength A^ this corresponds to a 
wavelength range A 4 (l + z min ) < A < A, 4 (l + z max ). 

4.1. Recovery of Lya optical depths 

We derive the Lya (1216 A) optical depth in each pixel 
from the flux F(X) using t L yQ (A) = — ln(F); pixels with 
Tya < are discarded. Pixels are considered "sat- 
urated" if F(X) < N a a\, where a\ is the rms noise 
amplitude at the pixel and N a is a parameter that we 
vary in our tests. For saturated pixels we can recover 
a good estimate of the true Lya optical depth using up 
to N\ w higher-order Lyman lines. We define r^ a = 

min{TLy„/LyaALya//LynA L yn}) where / Ly „ is the oscilla- 
tor strength of the nth order Lyman line and AL yn is its 
rest wavelength (Lyl= Lya, Ly2= Ly/3, etc.). We use all 
lines with 1 < n < N\ m that lie in the wavelength cov- 
erage of the spectrum and for which N a a n < F(\-L yn ) < 
1 — N a a n , where a n is the noise at ALyn- The pixel is 
discarded if none of the available higher orders satisfies 
this criterion. Taking the minimum optical depth min- 
imizes contamination by other lines, while the selection 
criterion picks out strong but unsaturated lines and mini- 
mizes the effects of noise. In particular, this criterion au- 
tomatically excludes regions in which the noise is greater 
than a > {2N a )- 1 . 

The recovery of saturated Lya optical depth is illus- 
trated in Figure [l]. For these and other tests we have 

http : //pitto . mib . inf n .it/ ~haardt/ref model . html 

plying the optical depth by a fixed factor) instead of the amplitude of 
als is included. 
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generated fiducial simulated spectra as described in § 
with z qso = 3.5 (top row) and z qso = 2.5 (bottom row 
wavelength coverage of 3000-7500 A, 0.04 A pixels and 
6.6 kms -1 detector resolution. A constant metallicity of 
Z = 0.01 Zq is assumed for the gas. Noise is added to each 
spectrum to give a signal-to-noise ratio S/N as described 
in § 0; for most tests we quote results with both S/N = 25 
and S/N — 100, to bracket plausible, high quality obser- 
vations. 

All panels plot the recovered Lya optical depth r rcc vs. 
the 'true' (i.e. free of noise, instrumental broadening, and 
metal-contamination) value Tt ruo . The grey lines indicate 
the 25th(75th) percentiles of r roc . 4 For the left two pan- 
els of the top and bottom rows of Fig. II] we use A CT = 3, 
iVh = 10 and S/N = 100 or S/N = 25. The optical depth 
can be recovered quite reliably for T\^ ya up to several hun- 
dred; the errors are only slightly larger for S/N — 25 (sec- 
ond panel of top and bottom rows). The main error is a 
slight over-estimation of optical depth for saturated pixels, 
due to contamination (as shown in the bottom four panels, 
the effect is much smaller at z = 2.5 where the Lya for- 
est is less crowded). Note that for fixed spectral coverage, 
the number of higher-order lines available for reconstruct- 
ing Lya PODs depends strongly on redshift: for coverage 
from A > 3000 A, all higher order lines are available only 
for gas at z > 2.29, and only two higher-order lines are 
available for gas at z m 2.1. The middle row shows trials 
with the same parameters as the top-left panel, but with 
Aho = 1, 2, 4, or 8; this illustrates both the effectiveness of 
the recovery and the importance of good spectral coverage 
of the higher-order lines. 

In the third panels from the left (in the top and bot- 
tom rows) we show N a = 1 (with A ho = 10, S/N = 25); 
lower N a gives less accurate recovery (the 'bars' extend- 
ing to rtrue > T roc are caused by noise), though less pixels 
are discarded as having no reliable optical depth estimate. 
Even for N a — 3 the number of discarded pixels is quite 
small (none at z = 3.5 and ~ 10% at z = 2.5), so we set 
N a = 3 in all subsequent trials unless otherwise noted. 5 

The rightmost panels (of the top and bottom rows) show 
the effects of errors in the continuum fit, estimated using 
the following procedure. We divide the spectrum into 90 A 
bins with centers Xk, and find for each bin the median flux 
fk ■ We then interpolate a spline across the fk , and flag all 
pixels that are at least N^a below the interpolation. The 
medians are taken again using the unflagged pixels and the 
procedure is repeated until the fit converges. The fluxes 
and errors are then scaled by the fitted continuum. The 
accuracy of the fit depends upon N£ a; fitting in highly- 
absorbed regions is more accurate with low A£ (because 
more low-flux pixels are rejected), whereas relatively un- 
absorbed regions are better fit with high A£ f . We find that 
using A^' = 1 when analyzing C IV and using A£ =0.5 
for O VI gives an acceptable continuum fitting error in 
the metal line regions (and either value gives unimpor- 
tant errors in the Lya region). Comparison between the 
original and scaled spectra then shows that for z qso = 3.5 



the induced errors are ~ +0.23 — 1-1.3% in the Lya re- 
gion, ~ +1.4% in the O VI region and ~ +0.06% in the 
C IV region; for z qso = 2.5 the errors are ~ —0.4 - +0.4%, 
~ 0.1%, and ~ —0.1% respectively. As shown in the figure, 
this continuum error does not have any significant effect 
on the recovered Lya optical depths; it will however, be 
more important in recovering the metal-line optical depths 
as discussed below. 

4.2. Recovery of metal line optical depths 

Having obtained for each pixel a good esti mate of the 
Lya optical depth (which, as shown below in § 5.1 can be 
converted into a density of the absorbing gas), we now turn 
to the metal-line optical depths for various species, which 
we will use to recover information about the metallicity of 
the absorbing gas. 

For each pixel with z m j n < z < z max we find the wave- 
length corresponding to absorption by some metal species, 
A = Aj(l + z), where A^ is the rest wavelength of the metal 
line transition; in this study we investigate the recovery of 
the species O VI, C IV, N V and Si IV, concentrating on 
the first two. The species we consider all appear in dou- 
blets (labeled here with subscripts T and 'II'), with the 
stronger transition (T) twice as strong as the weaker (i.e., 
fi,iK,i — 2fi,iiK,n where / is the oscillator strength). 
We then compute a metal line optical depth associated 
with each pixel of Tj(z) = min(Tij(z), 2nji{z)), i.e., the 
minimum of the optical depth of the stronger and twice 
the optical depth of the weaker component of the dou- 
blet. 6 Taking the minimum decreases contamination ef- 
fects, which are particularly important for N V and O VI. 
Different from previous work, we do not discard pixels for 
which neither component is well-detected. Instead, we as- 
sign an optical depth r^(z) = — \n(N a a\ i ,) (if the line is 
saturated) or Ti{z) = r m i n = 10~ 5 (if the optical depth is 
negative). Because there is little scatter in metal POD at 
high Lya POD, the former prescription causes tz(thi) to 
flatten smoothly to the highest detectable value. The lat- 
ter prescription allows pixels with very low optical depth 
to correctly influence the median POD (the number of 
well-detected versus poorly-detected pixels contains infor- 
mation about the metallicity) as long as T m - m is smaller 
than the median in any bin, which holds for T mul < 10~ 5 . 

Given z m - m > zp, the four chemical species we consider 
suffer very different levels of contamination from hydrogen 
lines. Carbon IV is contaminated only by other metals and 
by its own doublet, and Si IV is similarly uncontaminated 
through most of its available wavelength range. On the 
other hand, O VI is contaminated by several Lyman lines 
including Lya, and the wavelength range for N V overlaps 
almost completely with that of Lya (although there is a 
tiny segment with Sz ~ 0.1 for which N V lies red- ward of 
the QSO's Lya emission line and is thus relatively uncon- 
taminated) . The importance of contamination can be seen 
in Figs. H and ||, which show, for z qso = 3.5 and z qso = 2.5 
respectively, examples of a small spectrum segment in Lya 
region (top panel) and the corresponding C IV (middle 

4 In each bin the 25th(75th) percentile represents the median of pixels that are below(above) a curve linearly interpolating the bin medians. 
They should not be interpreted as errors on the medians. 

5 For a realistic noise array in which noise increases in the far blue, somewhat more pixels are discarded because coverage of the higher-order 
lines is lost. Particularly at z < 3, this might call for a lower value of JV CT to maintain good statistics at large r values. 

6 Self-contamination can be corrected (see below) in regions (such as C IV) where it is more important than contamination by other lines, 
and in this case the minimum should not be taken. 
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Fig. 1. — Tests of the Lya optical depth recovery for saturated pixels, using higher-order transitions 
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All of the saturated pixels and a random 10% of the rest are plotted, and the grey lines indicate 25th and 75th percentiles of 
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of noise, cont. 
described in § 

recovered optical depth in each bin of true optical depth (see text). For this test we assume z qso = 3.5 (top and middle panels) or z qso = 2.5 
(bottom panels). The left two panels of the top and bottom rows show results for 10 higher order transitions, a threshold N a = 3, and 
S/N = 100 or S/N = 25 in the simulated spectrum. The third panels from the left on the top and bottom rows use N a = 1. The rightmost 
panels on the top and bottom rows test the sensitivity to errors in the continuum (see text for details.) The middle row uses Zqso — 3.5, 
N a = 3, and S/N = 100, but its four panels show N^ a = 1,2,4, or 8. We use a uniform metallicity Z = 0.01 Zq, 0.04 A pixels, 6.6kms -1 
detector resolution, and a spectrum extending from 3000-7500 A, and show pixels redward of Ly/3 absorption. 
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panel) and O VI regions (bottom panel). All three pan- 
els are plotted both with and without contamination by 
Lyman series and metal lines. 

A significant fraction of the contamination in the OVI 
region is due to higher-order Lyman transitions, and this 
contamination can be partially removed using the (accu- 
rate) estimate of the Lya optical depth that has been com 
puted using those higher-order lines as described in § 4.1 



To do this, for each pixel with wavelength A in the O VI re- 
gion, we cycle through N corc higher-order lines. For each 
line i = 1, . . . , N corr we check the recovered Lya optical 
depth TL yQ at wavelength A' = (AL y a/AL y i)A, and subtract 
the optical depth Ty. a (A') (fhyi^hyi/ fhya^hya)- As can be 
seen in Figs, g and ||, this correction is quite helpful. For a 
given QSO redshift, the importance of this correction in- 
creases (as does H I contamination itself) with decreasing 
redshift of the absorbing gas, because progressively more 
higher-order Lyman lines contribute to the contamination. 

Figure [i] shows how well the true O VI optical depths 
can be recovered, using a simulated spectrum. Each panel 
plots the median binned recovered O VI POD versus the 
true one (i.e. free of noise, contamination, and instrumen- 



tal broadening). Top panels are for 



3.5 while bot- 



tom panels show z qso = 2.5. The left panels show the 
median and 25th and 75th percentiles in recovered rovi 
for our fiducial model (note that these percentiles are not 
errors in the median) . There is a clear correlation between 
the recovered and true O VI optical depths and the recov- 
ery works particularly well at z qso = 2.5 where there is far 
less contamination. We have assumed a uniform metallic- 
ity; the correlation signal would be steeper if the metallic- 
ity were to increase with density. The panels on the right 
show the median curves with more noise (S/N = 25), with- 
out correcting for contamination from higher-order Lyman 
lines, or with errors in the continuum fit included. Note 
that in both cases the continuum fitting errors work in the 
opposite sense from most others, decreasing t he O VI opti- 



cal depths; this effect is discussed further in § |5.3| . Overall, 
recovery of O VI PODs is imperfect, even for fairly large 
r ovi ! but in § ^| we shall discuss how the simulations can 
be used to correct for this imperfect recovery. 

Self-contamination can be corrected for if it dominates 
other contamination, and we implement this correction for 
C IV. For each pixel at wavelength A in the C IV region 
we check the optical depth (computed without taking the 
minimum of the doublet) at the wavelength AA//A// where 
Xj and A// > A/ are the rest wavelengths of the C IV 
doublet. Half this optical depth (the theoretical strength 
of the second doublet component) is subtracted from the 
initial POD estimate, and the process is iterated (and con- 
verges after about 5 iterations). 

The C IV self-contamination correction works very well, 
but can enhance the effect of other strong contaminating 
lines that may be present in the C IV region of observed 
spectra (e.g., Mg II), by erroneously subtracting a doublet 
component as if they were C IV absorption. If such lines 
are not removed by hand, it is therefore important to re- 
move these contaminants; this can be done automatically 
by discarding all pixels in the C IV region with optical 



depth r(A) satisfying 

exp(-r) + 3cr < exp[-2r(AA /7 /A/) - r(AA//A//)/2], 

where a is the r.m.s. noise of the pixels at A, AA//A//, and 
AA///A/. Such pixels have too much absorption to be the 
sum of any doublet (with a primary component of twice 
the stength at AA/ /A// and a primary component (with a 
doublet of half the strenth at AA///A/), and are very likely 
contaminated. The self-contamination correction should 
then be done. We have tested this procedure by setting 
the flux F — in a random fraction / of the pixels in 
the C IV region. For / = the contaminant removal does 
not change the recovered tcivj & n d with both contami- 
nating lines and self-contamination removed as described, 
the recovery is not significantly affected for / <J 0.2 (if 
no contaminant removal is done, but the minimum of the 
doublet components is taken, then the recovery is still only 
significantly affected for / > 0.1) This robustness of me- 
dian statistics is, of course, one of the advantages of the 
median POD method. 

Figure |^ shows tests of the CIV recovery. As shown 
in the left panel, the recovery is nearly perfect down to 
tciv ~ 10~ 3 at both redshifts. In the right panel, the dot- 
ted and dot-dashed curves show how the recovery becomes 
slightly less accurate with S/N = 25 or with a continuum 
fitting error included. Here, continuum fitting errors arise 
more due to noise than to significant absorption and do 
not tend to systematically reduce optical depth 7 ; instead 
they act as an additional source of random error that re- 
duces the accuracy of the rciv recovery. The dashed line 
shows the curve without the self-contamination correction 
(but with the minimum of the doublet components taken); 
this shows that for z qso = 3.5 and our assumed metallicity, 
self contamination is the limiting factor in the C IV POD 
recovery if not corrected for (as was the case in previous 
studies). 

5. TESTING AND INTERPRETING THE POD METHOD 

The previous section described the recovery of metal and 
Lya PODs from a simulated (or observed) spectrum, and 
tested the accuracy of this recovery process. In this sec- 
tion we use our simulations to draw connections between 
the recovered PODs and physical quantities such as gas 
density and metallicity. 



5.1. Interpreting HI PODs 

Previous studies - both numerical (e.g., Zhang et al. 
1998, Dave et al. 1999; Croft et al. 1997,1998) and an- 
alytical (Schaye 2001) - have shown that there is a fairly 
tight relation between the column density or optical depth 
of Lya absorption, and the density of the gas responsible 
for the absorption. Figure confirms that a tight relation 
between optical depth ana gas density holds on a pixel- 
by-pixel basis. All panels plot the median Lya optical 
depth weighted overdensity versus the recovered 7L yQ for 
Zqso = 3.5 (top panels) and z qso = 2.5 (bottom panels), us- 
ing Nh — 10 and N a — 3. The left panels show the median 
overdensity (dark lines) and the 25th and 75th percentiles 
(light lines) for S/N =100. The relation is very close to a 
power-law, except at very low optical depth, T^ ya <^ 10~ 2 , 
where it flattens off. Least absolute deviation fits (dashed 

This is only true if 7V£ f is carefully chosen; otherwise the continuum fitting can systematically increase or decrease the optical depth in a 
relatively empty region. 
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Lyoj region with metal contamination included (light line) and not included (dark line). The middle panel shows the corresponding region for 
C IV. The black line includes only the stronger doublet component, and is not contaminated by other metals (the black line is nearly invisible 
because the dark grey line falls almost exactly on top of it). The light line includes both doublet components and contamination (by silicon). 
The dark grey line line shows the (contaminated) spectrum corrected for C IV seif-contamination (see text). The lower panel's dark line is 
the O VI region, with all contaminants turned off. The light line shows the spectrum with H I and other metal absorption included, and the 
darker line shows the spectrum after removal (see text) of the higher-order Lyman lines. A uniform metallicity of Z = 0.01 Zq was used for 
all panels. 
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Fig. 4. — Tests of the recovery of the 'true' OVI pixel optical depth from the simulated spectra. The left panels show the median, 25th 
and 75th percentiles in recovered optical depth, versus the noise-free, perfect resolution and uncontaminatcd optical depth. The top panels 
assume Zq SO — 3.5 and [^mini^max] — [3.1,3.5]; the bottom assume Zq SO — 2.5 and [zmim^max ] = [2.1,2.5]. The left panels use S/N = 100. 
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Fig. 5. — As in previous figure, but for CIV. Here, the top panels assume z qso = 3.5 and [z m i n , z m ax] = [2.8,3.5]; the bottom panels assume 
Zqso = 2.5 and [z m i n ,Zmax] = [1.95,2.5]. The three variations shown are: S/N = 25 (dotted lines), with continuum fitting errors included 
(dot-dashed lines), and without correction for self-contamination (but with the doublet minimum taken) (dashed lines). 
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lines) yield power-law indices of 0.60 (z = 3.5) and 0.68 
(z = 2.5). Neglecting redshift space distortions, theory 
predicts a power-law relation of TLy Q <x j 2 - 76 -o.767^ -^gj-g 
the gas temperature obeys T = T a 5 1 ~ 1 (Hui & Gnedin 
1997) with 1^7^ 1.6. This predicts then 6 oc if yct 
with 0.5 < a <J 0.65, in agreement with our results. The 
existence and tightness of this relation shows that for rea- 
sonable observational errors the density of absorbing gas 
can be recovered to within ~ ±50% about half of the time 
for each pixel, and very reliably on a statistical basis. This 
will allow us to readily translate information gleaned us- 
ing pixel optical depths into information regarding gas 
of specific density (though with some uncertainty com- 
ing from uncertainty regarding the ionizing background). 
The right panels show medians only, for S/N = 25 (dot- 
ted lines), with perfect detector resolution and no noise 
(dashed lines), and with continuum fitting errors included 
(dot-dashed). These curves are all quite close except for 
optical depths lower than the errors (mainly due to noise). 

5.2. Interpreting metal PODs 

Given the tight relation between POD and gas density, 
the ratio of metal to Lya PODs is related to the ratio of 
metal ion to H I den sity (subject to the errors in the re- 
covery discussed in § |4.1j ) , which is in turn related (by an 
ionization correction) to the gas metallicity. One potential 
pitfall in this chain of inference is that both ionization ef- 
fects and the difference in the thermal widths of metal and 
hydrogen absorption may lead to slight offsets between the 
Lya and metal lines (see Ellison et al. 2000 for some dis- 
cussion). It is therefore important to test whether the ab- 
sorption by metals and hydrogen arises in the same gas on 
a pixel-by-pixel basis. Figure tests whether the O VI and 
C IV absorption in a given pixel is due to gas of the same 
density as the gas giving rise to the Lya absorption. The 
left and middle panels plot the gas overdensity weighted 
by the C IV and O VI optical depth respectively (denoted 
by <5ovi an d <$crv) against the Lya-optical depth weighted 
density 8m for z qso = 3.5 (top panels) and z qso = 2.5 (bot- 
tom panels). The correlation is good but not perfect. For 
high $hi the metal (and in particular the O VI) optical 
depth weighted densities are low compared with 5m . This 
effect is mainly due to the greater thermal width of the hy- 
drogen lines (tests show that the correlation becomes very 
tight if the detector resolution approaches the typical Lya 
line width or if the metals are given atomic weight unity 
so that their lines are as broad as Lya). This discrepancy 
will always be unimportant for low (log(5) <J 0.5) densi- 
ties, but could potentially smear the Z{$) relation for high 
densities. However, as we show in § ^, this does not appear 
to be a significant effect, at least not for the density range 
over which we can measure the metallicity (5z <^ 10). 

If the gas temperature is close to ~ 3 x 10 5 K, then col- 
lisional ionization becomes important for O VI, and O VI 
is therefore often regarded as a probe of hot gas. To reach 
such high temperatures, the gas needs to be shock-heated. 
At low redshift (z 1), gravitational accretion shocks 
may well heat a significant fraction of the baryons to tem- 
peratures T ;> 10 5 K, but the fraction of gravitationally 
heated gas is predicted to be much smaller at the redshifts 
of interest here (z J> 2) (e.g., Cen & Ostriker 1999). In- 
deed, we find that the fraction of shock-heated gas in our 



simulation is small enough for collisional ionization to be 
negligable. However, it should be noted that our simu- 
lations contain only gravitational shocks. If shocks from 
galactic winds keep a significant fraction of the high-z IGM 
in a hot phase, then our simulations would underestimate 
the fraction of collisionally ionized O VI. If this were the 
case then the inferred median O VI optical depth would 
of course still be correct. However, as mentioned in § SI 
the (incorrect) assumption of pure photoionization would 
then lead to inacurate conclusions concerning the oxygen 
abundance. 

5.3. Testing the method 

Having shown that the recovered Lya optical depth is 
related to the gas density and that the metal absorption 
arises in gas with roughly the same density as the gas re- 
sponsible for the hydrogen absorption, we will now inves- 
tigate the relation between the measured T r £ c {r^ a ) and 
the metallicity as a function of the gas density. In this 
section we will investigate the effects of variations in the 
recovery procedure and the metallicity Z(5) on the cor- 
relation between metal and Lya optical depth. Figures || 
and ^ show the basic results of the POD method applied to 
simulated spectra with uniform metallicity Z = 0.01 Zq. 
Figure || applies the method to O VI, using fiducial pa- 
rameters (iVho = 10, iV corr = 4, and N a — 3). The top 
panels use [z m [ n , z max — z qso ] = [3.1, 3.5] while the bottom 
panels use [z m - m , z max = z qso ] = [2.1, 2.5]. Each panel plots 
T ovi vs - T Lya' binned into 20 bins of log(r£y C Q ) in the range 
[—3, 3] . Dark lines indicate medians in each bin, while light 
lines (in left panels) signify 25th and 75th percentiles. The 
medians and percentiles are calculated separately for each 
of ten independent spectra, and the average of the medi- 
ans and percentiles are plotted; the error bars about the 
medians in the left panels represent the standard deviation 
of the medians about the mean median in each bin. 

The correlation between t vi and TL ya r<=<=, which is 
clearly evident in the z qso = 2.5 case for r£y C a ;> 0.03, in- 
dicates the detection of O VI in gas of overdensity 5 ;> 0.2 
(using Fig. ^|); at z qso = 3.5 the correlation is significant 
at T£y C Q f> 0.1, or 5 ;> 0.5. The dashed horizontal line in 
each panel on the left represents the median T e yj, irre- 
spective of t£° c q . If no O VI were present, all data points 
would (and tests show do) lie near this horizontal line. 
Note that unless O VI contributes only negligibly to the 
mean optical depth in the metal region, this line is not a 
"reference level" (or detection limit) that gives the contri- 
bution of noise and other errors in the absence of detected 
metal absorption; the latter would correspond to the con- 
stant Tovi to which the curve asymptotes as r£° c a — > 0. 
One might hope to subtract such a "reference level" (as in 
Cowie & Songaila 1998) to see more clearly where metals 
are detected, but this requires a very accurate estimate 
of T e VI , and will generally lead to negative optical depths 
and large errors at low T^ ycc . Thus, for practical purposes 
we prefer not to subtract the reference level, and instead 
recover metallicity information as described below in § [| 

The right panels show the effects of different assump- 
tions or parameters in the analysis. The dashed lines, 
where no correction for higher-order Lyman lines in the 
O VI region was performed, show that this subtraction is 
quite helpful, and crucial if O VI is to be detected at z > 3 
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Fig. 6. — Lya optical depth weighted overdensity of gas responsible for the absorption in each pixel, versus the recovered Lyo POD in 
that pixel, binned in r£ ec . Left panels plot the median and 25th and 75th percentiles of gas overdensities corresponding to a given pixel 

optical depth; top panels show [z m i n , z max , z qso ] = [2.8,3.5,3.5] while bottom panels show [z m i n , 2 m ax, Zqso] = [1.95,2.5,2.5]. Left panels show 
the fiducial S/N = 100 model. A least-absolute-deviation power law fit is also included and gives the relations <5 = 1.2t®®® (z = 3.5) and 

<5 = 2.5r°^| (z = 2.5). Right panels (showing only medians) show trials with S/N = 25, with continuum errors included, and with no noise 
and infinite detector resolution. The optical depths are recovered using up to 10 higher order lines, and N a = 3. 
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Fig. 7. — Comparison of overdensity weighted by TL ya versus that weighted by tovi (left panels) or tqiv (middle panels), assuming 
Zqso = 3.5 (top panels) or z qso = 2.5 (bottom panels). The right panels plot the two metal optical depth weighted densities against each 
other. Metallicity is uniform at 1% solar, and there is no noise in the spectrum. 
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(this correction has not been employed in previous stud- 
ies). Likewise, as shown by the dot-dashed line (where the 
weaker doublet component is not used), taking the minu- 
mum of the doublet improves the optical depth recovery 
considerably. Additional noise (S/N — 25, dotted lines) 
degrades the correlation at low r^ a for z qso = 2.5, but has 
little effect for z qso = 3.5. Probably the most dangerous 
effect is induced by a significant error in the continuum 
fit, as this appears to generate a false correlation by sys- 
tematically decreasing O VI optical depths at small Lya 
optical depths (especially for z qso = 3.5 for which contin- 
uum fitting errors are large; see 4.1). This, however, is 



not as problematic as it first appears, because this effect 
will only slightly enhance an already-existing tovi vs. TL ya 
trend; we have verified that it will not create one. 

Figure ^| shows the same quantities for C IV. Again, C IV 
is clearly detected at very low density for Z = 0.01 Zq\ in 
this case metals are detected at t£° c q ;> 0.06 (S ;> 0.4) for 

= 2.5. 



z qso = 3.5, and at r£° y c Q £ 0.03 (5 ;> 0.2) for ^ qs 
Here, as in Fig. [s] above, the continuum fitting error be- 
comes a significant random error source at very small tqiv 
for z qso = 2.5, but the recovery is limited by noise for 
S/N < 50. The dashed lines show that correction for 
self-contamination (also not done in previous studies) is 
important at low optical depths. 

The results of the POD analysis applied to Si IV and 
N V arc shown in Figure EoL The left panel shows Si IV, 
with the usual choices of z qso , and with z m - ln chosen so 
that only Si IV absorption redward of Lya is used. De- 
spite the lack of Lya contamination, silicon is detected 
only at T^ a ^> 3 (or 5 ;> 3). The recovery is limited pri- 
marily by C IV contamination. Nitrogen V, shown in the 
right panel, can be detected well down to t£™ ~ 1 for 



■^qso 



2.5 but is nearly undetectable for 



'qso 



3.5. The 



dotted line shows the analysis using z min = 2.44 so that 
only the tiny region of N V absorption redward of Lya 
is used. The correlation is strong, but the statistics are 
poor and this region is probably too close to the QSO to 
be of practical use. Thus, while useful information can be 
gleaned from the analysis of Si IV and N V PODs, they 
are not recoverable at as low Lya optical depths as C IV 
and O VI. 

The overall aim of the POD method is to recover from 
spectra the abundances of various metals in low-density 
gas. Figure |ll| illustrate its effectiveness. The left pan- 
els re-plot the medians and percentiles for O VI (top) and 
C IV (bottom) absorption from Figs. || and|| as solid lines. 
Dashed lines show the results for a model in which the 
median metallicity of the gas is still 0.01 Zq, but with 
1 dex of scatter: the simulation is divided into cubes 
0.6/i _1 comoving Mpc on a side, and each is assigned a 
random metallicity drawn from a uniform distribution in 
—3 < \og(Z/ Zq) < —1. Compared to the uniform metal- 
licity case, the medians are very similar (though slightly 
higher) and the scatter slightly larger. This indicates that 
the POD method can determine the overall median metal- 
licity fairly robustly (see also Ellison et al. 2000) , but can- 
not yield strong information about the scatter in the metal- 
licity. 

The right panels show trials in which the metallicity 
varies with gas density. Solid lines are as in the left panels, 
and show results for a uniform metallicity Z = 0.01 Zq. 



Dotted lines are trials with Z = 0.01 Zq for 8 > 5 and 
Z = for S < 5. Here the Tz curve flattens, for both 
metals, at log(TLy C a ) « — 0.3 (or S 4 — 6), con- 
firming that the correlation vanishes at HI optical depths 
roughly corresponding to densities at which there are no 
metals. The difference between the solid and dotted lines 
at 0.2 <; log(r£y C Q ) ;> 0.7 does, however, indicate that the 
presence of metals at S < 5 affects metal PODs corre- 
sponding to HI PODs with log(rLy a ) > 0. The dashed 
line shows a trial with Z — 0.001 Zq. For carbon the 
curve is quite similar to that for Z = 0.01 Zq, scaled down 
by ~ 1 dex; for O VI the signal is quite washed out even 
for larger TL yQ . Finally, the dot-dashed line is for metal- 
licity Z — 0.0015 Zq. For C IV this curve roughly matches 
the Z = 0.01 Z Q curve at log(T LyQ ) w 0.8 (S « 13), 
and matches the Z — 0.001 Zq curve for log(rL yQ ) ~ — 1 
(6 ~ 1), as one would hope (though in the latter case the 
values are determined by the noise). The trend of increas- 
ing metallicity with gas density is apparent in the steep 
slope of the dot-dashed line in the C IV plot, but the dif- 
ference in slope between this model and the constant 1% 
metallicity model with a cutoff is not very apparent for 
O VI. 

The results of this section show that the interpretation 
of the median pixel optical depth is reasonably straight- 
forward as long as noise and other systematic errors are 
small (as for C IV) but becomes somewhat more ambigu- 
ous when contamination or self-contamination are impor- 
tant (as for O VI). In the next section we describe how 
simulated spectra, which map "true" to "recovered" opti- 
cal depths, can be used to directly invert the pixel optical 
depths into C IV/H I and O VI/H I ratios as a function of 
gas density. 

6. RECOVERING Z(5) 

In the preceding sections, we described, tested, and im- 
proved the basic POD recovery method that has been used 
previously in analyses of absorption spectra. Despite our 
improvements, these recovered PODs are still somewhat 
affected by contamination, noise, etc. With our simulated 
spectra it is, however, possible to extend the basic method 
by using the relation - derived from the simulations - be- 
tween true and recovered PODs, to recover a good estimate 
of the true metal PODs from observations. Furthermore, if 
we use the simulations to determine the relations between 
density and Lya optical depth, and between density and 
temperature, and if we assume a fixed ionizing background 
radiation, then we can correct for the ionization of both the 
metal and hydrogen and convert the ratio Tz /t Lya (r Lya) 
into a metallicity as a function of density. In this section 
we describe and test this procedure. 

6.1. Correcting the optical depth using simulations 

Figures |^ and [| show examples of recovered vs. true 
pixel optical depths in O VI and C IV, for particular sets 
of parameters describing the simulated spectrum, and for 
a metallicity Z(S, z) = const = 0.01 Zq of the gas as 
a function of overdensity and redshift. Given the rela- 
tion Tjj oc (T| ruo ) for a spectrum simulated to have the same 
properties as an observed spectrum, the binned recovered 
metal PODs (such as those shown in figs. |^-[TT|) can be "in- 
verted" to yield "true" PODs, significantly removing the 
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Fig. 8. — Recovered OVI optical depth vs. recovered Lya optical depth for z qBO = 3.5 (top panels) and z qBO = 2.5 (bottom panels). Left 
panels assume a uniform mctallicity Z — 0.01 Zq , ^ m in — 3.1, ^max — -Sqso % No- = 3, and S/N = 100. The median and 25th and 75th 
percentiles are averaged over ten independent spectra, and the error bars represent the standard deviations of the medians about the mean 
median. The horizontal dashed lines indicate the median rovi; irrespective of r^ ya . Right panels show mean medians for the fiducial model 
(solid lines) and various other models: with S/N = 25 (dotted lines), without correction for higher-order Lyman lines (dashed lines), without 
taking the minimum of the doublet, (dot-dashed line), or with continuum fitting error included (triplc-dot-dashed lines). 
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Fig. 9. — Recovered CIV optical depth vs. recovered Lya optical depth for z qso = 3.5 (top panels) and z qso = 2.5 (bottom panels). Left 
panels assume uniform metallicity Z = 0.01 Zq, z m ; n = 3.1, z max = z qso , A^c = 3, and S/N = 100. The median and 25th and 75th percentiles 
are averaged over ten independent spectra, and the error bars represent the standard deviations of the medians about the mean median. The 
horizontal dashed lines indicate the median tqiVj irrespective of iXya- Right panels show mean medians for the fiducial model (solid lines) 
and various other models: with S/N = 25 (dotted lines), without self-contamination correction (dashed lines), or with continuum fitting 
errors included (dot-dashed lines). 
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effects of noise, thermal broadening, contamination by H I, 
and continuum fitting errors. The remaining effects of self- 
contamination and contamination by other metal lines, if 
important, can also be corrected for provided the metallic- 
ity used in the simulated spectra (which are in turn used 
to perform the inversion) is similar to the metallicity in 
the observations. 

Figure |l2| exhibits tests of the importance of metal line 
contamination for the recovery of true optical depths for 
O VI (left panels) and C IV (right panels), for z qso = 3.5 
(top panels) and z qso = 2.5 (bottom panels). Curves of 
r z cc ( T z ue ) are plotted for several different choices of Z{S). 
With the exception of O VI at z = 2.5, the curves are vir- 
tually identical in their overlap region, indicating that the 
deviation from perfect recovery is dominated by factors 
other than metal line contamination (O VI at z = 2.5 has 
a significant level of self-contamination; recall that self- 
contamination was removed for C IV). If the curves are in- 
dependent of metallicity within the range of c recovered 
from the analyzed spectrum, the inversion can be done in 
one step; if the assumed metallicity affects the r^ ec (r^ rue ) 
obtained, the inversion can still be done using an iterative 
approach, as discussed below. 

The basic inversion procedure is demonstrated in 
Fig. Ol which illustrates the process applied to C IV in 
a Zqso — 2.5 spectrum with fiducial parameters. Given 
T zk ( T hya k) m kins k (panel 1) and the monotonic (hence 
invertiblc) curve t^ c (t^ uc ) fit though the binned pix- 
els (panel 2), we are able to compute 7j^X( r Lya k) = 
r | rue t T Ffc ( r Lya fc)] (hereafter we drop the index k of the 
bin). The t™ c (t|J' uc ) relation, which we compute by aver- 
aging over 10 spectra, is indeed monotonic but increasingly 
flat at low ri rue . We therefore determine a value 

Zj ^,111111 

(and T^'min ) (shown in light, dotted lines in panels 2 and 



3) below which the curve is effectively flat, and determine 
an "error" cr m i n on this value by computing the standard 
deviation of the median optical depth across the ten real- 
izations for all of the pixels with "r| rue < r^'^ in . The curve 
also varies from one realization of the spectrum to the next; 
these variations are quantified in each bin by the standard 
deviation of the 10 averaged realizations, and two more 
monotonic functions T^'^ir 1 ^ 6 ) can be constructed us- 
ing the ±lcr values. The median, ±ler, and ±2ct values of 
r!j ec are then inverted into 7j? v : the median values are in- 
verted using the mean r^ ec (t^' uc ) , while the — 1(7 and — 2er 
values are inverted using the r™ c ' + (rf uc ) curve and the 
+ 1(7 and +2<t values are inverted using the T^ ec '~ (r^ ruc ) 
curve. As for the "flat" part of the recovered spectrum, 
we set t|? v = when r^ c < T^^ in , and set the ±ler and 
±2ct values to zero when they fail to exceed r^ e ^ Jn =p oWn- 
Projecting the errors this way ensures that: A) we have 
included a conservative estimate of the error induced by 
variance in the true vs. recovered POD relation from one 
realization of the spectrum to the next, and B) the flat- 
tening of t^ >c (t^' uc ) correctly manifests as large negative 
errors (effectively upper limits) in r!? v for low values of 

rcc 
T hya- 

The result of this procedure is shown in panel 3. The 
plotted points are our best estimates of the true, uncon- 
taminated tz versus TLya ■ These will be reliable to the 
degree that errors in the recovered PODs are dominated 
by noise or H I contamination, which should be correctly 
modeled by the simulations. If the errors are dominated 
by continuum fitting errors, the inversion procedure should 
correct for this if the same continuum fitting procedure is 
used (with equal effectiveness) on the simulated and real 
spectra. If (as for O VI at \ow-z) self-contamination is 
important, the inverted PODs will be accurate only if the 



19 




CIV, 



,= 2.5, S/N=100 




Median, u = dex 

Median, (7=1 dex — — — 

25, 75%, (7 = dex 

25, 75%, (7=1 dex 



-3 -2 -1 

log(Recovered T Lya ) 



o 



-1 



-2 





i 





: Z = 0.01 Z 




/' ' 


: Z=0.01 Z s , 5 cut = 


= 5 




-_ Z=0.001 Z Q 




- 


: Z = 0. 0015 Z 


/ / 


i : 




/ /: 


/ - 




/■' 


/ : 




/ ■■' ' 


















-3 -2 -1 

log(Recovered T LyK ) 



Fig. 11. — Recovered metal optical depth vs. recovered Ly« optical depth for OVI (top panels) and CIV (bottom panels) for different 
mctallicity distributions. Left panels use a median metallicity Z = 0.01 Zq, but one trial with a uniform metallicity (solid lines), and one in 
which there is a 1 dex scatter in log(Z) on ~ 1 Mpc scales (dashed lines). The right panels show trials in which the metallicity is constant at 
Z = 0.01 Zq (solid) or Z = 0.001 Zq (dashed), a trial in which Z = const, for overdensity <5 > <5 cu t = 5 and zero otherwise (dotted), and a 
trial in which Z = 0.0015 (dot-dashed). All panels have 2 qso = 2.5, S/N = 100, N a = 3, N ho = 10 and iV corr = 4. 




Fig. 12. — Test of the recovery of the pixel optical depth from the simulated spectra, for different mctallicitics. Left panels show OVI, while 
right panels show CIV. Top (bottom) panels have z qBO = 3.5(2.5). Recovered vs. true POD is plotted for uniform metallicities of 10 -3 Zq, 
10~ 2 Zq, and Z = 0.0015 .Zq. The light solid line indicates perfect recovery. The flattening at high tqiv is due to saturation. 
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level of metal contamination in the spectra used in the 
inversion is similar to that in the real spectrum. 

6.2. Interpreting the optical depth ratios 

As mentioned in § ||, the simulations can be further used 
to interpret the PODs in terms of the metallicity Z(S, z) 
of the absorbing gas. 

First, the ratio of the metal to H I optical depth can be 
converted into the ratio of the metal ion to H I density, 
both as a function of the H I optical depth. For example, 
for C IV we can compute for each bin, 



nciv 
n H i 



(/A) 



HI 



'Ly 



« (/A) 



CIV 



Second, the tight relations between gas density and Lya 
optical depth, and between temperature T and density, 
give for each bin in Lya optical depth, a density and tem- 
perature of the corresponding absorbing gas, denoted re- 
spectively by Sh ya (Tl c y c a ) and Tl vq (t l ° c q ). These are shown 
(solid lines) in panels 4 and 5 for our example inversion. 
Given an assumed ionizing background, these can then 
be used to compute, for each bin, an ionization correc- 
tion n#(nHi) to yield a physical H I density of absorbing 
gas for each bin. Third and similarly, the binned, metal 
optical-depth- weighted density and temperature, Sz(ij^ a ) 
and Tz(rl^ a ) (panels 4 and 5, dotted lines) can be used to 
correct for the ionization of the metal, to yield an estimate 
of the density rii of the atomic species i (i — C, O, etc.). 8 

The ratio ni/nu then represents a true metallicity of 
species i for each bin, yielding (using Sz[tz° c ]) an esti- 
mate of Z(8z)- The metallicities in our sample inversion 
are shown in panel 6, with the dark dashed line repre- 
senting the true metallicity. The good agreement between 
the measured and true metallicities demonstrates the effec- 
tiveness of the procedure (even when the metallicity used 
to derive t^ c (t^' uc ) is different from the true metallicity) 
when the ionizing background is taken to be the same as 
used in the spectrum generation. 

The solid and (light) dashed lines indicate the changes 
in recovered metallicities induced, respectively, by an over- 
all scaling of ±50% in temperature or density used in the 
ionization correction; these are conservative estimates of 
how much the temperature and density in the simulations 
could differ from those of the true IGM. Such variations 
induce only small changes in the inferred metallicity com- 
pared with the statistical errors. 

More uncertain is the shape and spatial variation of the 
ionizing background, and the importance of collisional ion- 
ization (collisionally ionized gas would have very differ- 
ent ionization corrections than photoionized gas); here we 
have used the model of the UV-background from quasars 
and galaxies of Haardt and Madau (2001), and used sim- 
ulations without feedback from star formation (i.e., only 
gravitationally induced shocks are included) so that the 
fraction of collisionally ionized gas is negligable for the 
redshifts of interest here. The effects of changes in the 



assumed ionizing background radiation can be studied by 
varying the background used in the inversion, but we will 
not do this here. 

After the recovery of metallicities Z(6), the inversion 
can be checked or if necessary iterated, as follows. A new 
set of simulated spectra can be generated using the metal- 
licity found in the inversion, 9 and used to generate a new 
r z SC ( T z Ue ) (P ane l 2). Using this new relation the inver- 
sion can be repeated. If the results are consistent with the 
first iteration, then the recovery was not significantly af- 
fect by metal-line contamination. If the first two iterations 
differ, the procedure can be repeated until convergence is 
obtained. 

After convergence, the measured metallicity Z[S) can 
again be used to generate new simulated spectra. Com- 
parison of the pixel statistics (e.g., T v £ c (T™y a )) for these 
spectra and the original spectra then provides a strong 
consistency check on the recovered metallicity, the simu- 
lation assumptions, and the observed spectrum. 

The top two panels of Fig. [l4] show the results of 
the inversion procedure applied to C IV for z qso = 3.5 
(left panel), and for z qso = 2.5 with a true metallicity 
Z = 0.0015 Z G (right panel). The first shows that for 
Zqso — 3.5 the inversion is able to constrain the metal- 
licity of even underdense gas (we find that this holds for 
Z ;> 0.001 Zq); the second demonstrates that the inver- 
sion can cleanly recover a trend in Z{8). (Note that the 
^qso = 3.5 examples use only about half of the available 
rcdshift range; it is worth noting that here - and in general 
- this does not significantly increase the error bars.) The 
bottom two panels show the inversion applied to O VI, for 
z qso = 2.5 (left panel) and z qso = 3.5 (right panel). De- 
spite contamination by H I, O VI serves as a good tracer of 
metallicity for z qso = 2.5, even in slightly underdense gas. 
For z qso = 3.5 the inversion is not as accurate as for C IV, 
but can still yield a good constraint on the metallicity of 
overdense gas. 



7. DISCUSSION AND CONCLUSIONS 

Studies of the absorption spectra of QSOs at 2 <; z <J 4 
have revealed that the intergalactic medium is enriched 
with metals to the level of ~ 0.1 — 1% solar metallicity at 
overdensities S > 5. However, a more accurate estimate 
of the metallicity and its spatial variation, as well as the 
extent to which it varies with gas density and redshift, are 
currently unknown. This paper is the first in a series that 
systematically analyzes a sample of QSO absorption spec- 
tra in an effort to glean as much information as possible 
about the enrichment of low-density gas at high redshift. 

In the present paper, we have systematically described, 
tested, and discussed the "pixel optical depth" (POD) 
method of analyzing metal absorption in QSO spectra. 
This method, while still less widely used than direct line- 
fitting techniques, is better able to extract information 
from low-density regions in which individual metal lines 
are difficult or impossible to detect. The POD technique 

8 We find that using Lyct-weighted quantities to correct for H ionization and metal-weighted quantities to correct for metal ionization works 
slightly better than using either set of quantities in both corrections. 

9 In practice this may be done by performing a power-law or geometric fit to the recovered Z(&). If self-contamination is important, the 
inversion is inaccurate if the metallicity assumed in performing the inversion is very different from the true metallicity. This effect appears 
chiefly for tz below the dashed line in panel 1 (since, as mentioned in § |5 ■ 3| this line indicates the level below which the metal contributes 
significantly to the absorption). Thus, in the first iteration we fit Z(8) only using points with Tz above this dashed line. 
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Fig. 13. — Test of the recovery of gas metallicity versus density. Here, z qBO = 2.5 and we recover carbon. The analyzed spectrum has 
Z = 0.01 Zq, whereas the spectra used in performing the inversion use Z = 0.0015 Zq. Panel 1: r^?^ versus (diamonds). 1- and 2-cr 

errors determined using bootstrap resampling are shown. Panel 2: The median binned 'true' vs. 'recovered' metal POD relation. Errors are 
computed by taking the standard deviations among 10 independent spectra. The horizontal (vertical) light dotted line indicates the recovered 
(true) tciv value below which the relation is taken to be flat. Panel 3: 'Inverted' metal PODs using panels 1 and 2, as described in text. 
Dotted line is as in panel 2. The dashed line re-plots the relation of panel 1 for comparison. Panel 4: Lyo (solid line) and metal (dashed 
line) optical depth weighted ovcrdensity of absorbing gas vs. recovered Lyc* POD. Errors are computed as in panel 2. Panel 5: As panel 3 
but for temperature. Panel 6: Recovered carbon abundance in solar units versus metal optical depth weighted density, given an ionization 
correction using density and temperature from panels 3 and 4, and the same ionizing background as was used in the generation of the spectra. 
The true metallicity is plotted (dashed line) for comparison. Solid and light dashed lines indicate, respectively, changes induced by an overall 
±50% scaling in temperature or density in the ionization correction. 
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Fig. 14. — As panel 6 of the previous figure, but for different elements or metallicitics. Top left: CIV inversion for z qso = 3.5, with 
Zmin = 3.3. Top right: CIV inversion for z qso = 2.5, but with Z = 0.0015 Zq in the analyzed spectrum, and Z = 0.01 Zq in the spectra 
used for the inversion. Bottom left: OVI inversion for z qso = 2.5. Bottom right: OVI inversion for z qso = 3.5, with z m ; n = 3.3. 
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was developed and previously employed in QSO absorp- 
tion studies of C IV (Cowie & Songaila 1998; Ellison et 
al. 2000) and O VI (Schaye et al. 2000a), but has suffered 
from difficulties in the interpretation of the results. Here 
we overcome this difficulty by testing the method on spec- 
tra generated from realistic cosmological hydrodynamical 
simulations that are able to reproduce the observed statis- 
tics of the Lya forest. 10 This has allowed us to both assess 
the effectiveness of the method, and to refine and calibrate 
it. The major improvement to the method developed and 
tested here is the removal of higher-order Lyman lines from 
the O VI region, and the correction of self-contamination 
of C IV by its own doublet. These innovations significantly 
improve the accuracy of the O VI and the C IV recovery. 
In the Appendix we gave a compact but complete recipe 
for implementing the method. 

The most general result of our tests is that the POD 
technique is very effective at recovering information about 
the metallicity of low-density gas, even where (as for O VI) 
the metal lines are severely contaminated. More specifi- 
cally, our tests reveal the following: 

1. Lya optical depths of up to several hundred can be 
reliably determined using higher order Lyman lines, 
and these optical depths in turn give an accurate 
estimate of the density of gas giving rise to the ab- 
sorption. 

2. With high-quality spectra (signal-to- noise ;> 50) me- 
dian C IV optical depths binned in Lya optical depth 
can be accurately recovered down to tciv ~ 10~ 3 for 
both z qso = 3.5 and z qso = 2.5. Previous analyses 
were probably limited by self-contamination; with 
our correction for this effect, the recovery is limited 
by noise and by errors in the continuum fitting. 

3. The recovery of median O VI PODs is limited pri- 
marily by contamination due to H I lines (and sec- 
ondarily by continuum fitting errors). A significant 
fraction of this contamination can be removed, al- 
lowing a fairly accurate recovery of O VI PODs with 
Tovi 10~ 2 for z qso = 2.5, and a useful recovery of 
median PODs for z qso = 3.5. 

4. The Si IV (limited primarily by C IV contamina- 
tion) and N V (limited primarily by Lya contamina- 
tion) PODs can be usefully recovered for Lya optical 
depths (r LyQ ^ 1). 

5. Using the POD method both C IV and O VI should 
be detectable in realistic spectra of z ~ 3 QSOs even 
in gas of the mean cosmic density, if the metallic- 
ity is ~ 1O~ 3 Z , and at even lower density if the 
metallicity is higher. 



6. When applied to simulated spectra with different 
metallicities as a function of gas density, the method 
is able to distinguish models in which the metallicity 
is constant from those in which the metallicity de- 
clines significantly with gas density. The results are 
not very sensitive to spatial variations in metallicity 
at a given gas density. 

7. Most important for accurate recovery of metallicity 
information are large wavelength coverage (so that 
a significant number of higher order Lyman lines 
can be used), accurate continuum fitting, and high 
signal-to- noise. However very high S /N is only useful 
for C IV recovery; obtaining a S/N ^> 50 probably 
does not significantly improve the amount of infor- 
mation that can be recovered regarding O VI. 

8. Due to different levels of thermal broadening in met- 
als vs. hydrogen, the inferred density/temperature 
of the gas is, somewhat different when weighted by 
metal line optical depth than when weighted by Lya 
optical depth, especially at relatively high gas den- 
sities. However, this uncertainty does not appear to 
significantly affect the derived metallicities. 

9. Using spectra generated from cosmological simula- 
tions, the POD method may be extended, as de- 
scribed in § ^, to directly and accurately recover the 
metallicity of intergalactic gas vs. its density. The 
largest uncertainty in this "inversion" is the spec- 
tral shape (and possible spatial variations) of the 
ionizing background. For a given background, the 
abundances of C and O can be measured with errors 
of at most a factor of a few over at least an order 
of magnitude in density, using a single high-quality 
spectrum. 

The POD technique, as refined and calibrated in our 
study should, when applied to QSO spectra of currently- 
available quality, yield clear information about the metal- 
licity of even underdense gas. This information will place 
strong contraints on models of the enrichment of the IGM, 
and consequently on feedback and galaxy formation. 
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APPENDIX 

SUMMARY OF THE METHOD 

The POD method as we have described and tested above can be implemented as follows, given an absorption spectrum 
of a QSO at redshift z qso . 

1. Choose redshift range. Determine zp = (1 + z qS o)(AL y/ 3/AL yQ ) — 1, and choose a redshift range zp < z m i n < 
z max < z qso . For a given transition i with rest wavelength Ai , pixels with wavelength Aj(l + z m j n ) < A < Aj(l + z max ) 
can then be analyzed without Ly/3 contamination of the Lya region. 

2. Derive Lya PODs. Take TL ya (z) = — ln(F) where F(X) is the observed normalized flux at A = AL yQ (l + z); mark 
as "discarded" pixels with TL yQ < 0. 

If F(X) < 3a\ (where <j\ is the normalized noise array), the pixel is saturated. For these pixels, find the minimum 
of the Lya optical depth corresponding to the optical depth seen in the available higher-order lines: define 

T hy°a = mm { T LynS'Lya/SLyn} , 

where 5L y n = fhyn^hyn, fhyn is the oscillator strength of the nth order Lyman line and AL yn is its rest wavelength. 
Compute this, using as many higher-order lines as fall within the spectrum's coverage, use only pixels that are not 
"poorly detected", i.e., use those for which 3cr„ < F < 1 — 3<r n (where a n is the noise at AL yn )- If no higher-order 
lines are usable, mark the pixel as discarded. 

3. Derive metal PODs. Consider the spectral region covering all multiplet components of the metal in question for 
redshift z m j n < z < z max . Then find the optical depth in all multiplet components, that is, compute Tz,k{z) by using 
the flux at wavelengths Az,fc(l + z ) where \z,k is the rest wavelength for the fcth multiplet component of the metal 
species (relevant lines appear in doublets with k = 1,2). If the flux F < io\ z fc (i+z), the line is effectively saturated; 
in that case set Tz{z) = — ln(3er Az k (\+ z ))- 

A. Correct the metal PODs for higher-order Lyman contamination. If the metal lines lie in a region 
shared with higher-order Lyman lines (as is the case for O VI), subtract this contamination using an interpolation 
of the r£y C a values derived in step 2, as follows (otherwise skip to step C): 

TZ,k{z) ■= T Z ,k(z) - ^^ T L y a( AA lWA Ly n), 

where A = Az i fe(l + z), and the sum is over the higher-order lines with wavelength covered by the spectrum. Do this 
for each multiplet component k. Skip now to step C. 

or B. Correct the metal PODs for contamination and self-contamination. For lines without Lyman-series 
contamination, and in which self-contamination is dominant (such as C IV), self-contamination may be corrected 
using the following procedure. 

First correct for contamination by other strong lines: find all pixels with 

exp[-r z (A)] + 3ct < cxp[- (52/51 )r(X s ) - (gi/ 92)^X4)], 

where g k = fz,k^z,k, A = X z ,i{l + z), X d = X z ,2^ + z), X s = {X z ,i/X z ,2)X, and a = (a 2 Xg + a\ d +cr^) 1 2 is the 
quadrature sum of the noises at A s , A^, and A. Since the optical depth at A should be the sum of a primary component 
(with a doublet of strength 52/51 at A^) and a doublet (with a primary component of strength 51/52 at A s ), the 
selected pixels are then very likely contaminated; discard them so that they are used neither in the analysis nor in 
the self-contamination correction. 

Then, correct for self- contamination using the following algorithm: 

1. Begin with the optical depths of the primary component: 

t| cc (A)-t(A), A = A z ,i (! + *)■ 

2. Subtract the doublets corresponding to these optical depths: 

r| oc (A) := r(A) - (52/51)^^,1^,2). 

(The values of t t £ c at the specified wavelengths must be interpolated from those available; for wavelengths lying 
outside of the available range, take optical depths from the original full spectrum.) 

3. Repeat the second step until convergence occurs (approximately four more iterations). The resulting t v £ c (z) is 
well-corrected for self-contamination. 

C. Take minimum of doublet If self-contamination is not being corrected (i.e. step B was skipped), take 
t~z(z) — min{5 ; 7 Tz,k(z)} using only those components k > 1 for which Tz.k > 0. As above, 51 = 1, gk/gi = 
fz,k^z,k/.fz,i^zs, and fz.k is the oscillator strength (for O VI, C IV, N V and Si IV, 51 = 1, 52 = 2). If tz{z) < 0, 
set Tz(z) to some fixed value (e.g., 10~ 5 ) much smaller than the noise. 



27 



4. Bin the Lya and metal PODs. Create bins of —3 <; log(r£y C Q ) <J 2 and compute the median T^ ya and n in each 

bin. Errors can be computed by bootstrap resampling: divide the analyzed region into 5 — 10A chunks (the chunk 
size should be greater than the line widths), then compute a large number of realizations in which enough chunks 
are randomly selected (each chunk may be selected multiple times) and concatenated to form a spectrum with a 
number of pixels equal to the original. Bin each realization, and take as the error in each bin the standard deviation 
from the mean median. 



